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Objectives: To evaluate the association between fracture risk and levothyroxine use in elderly women with hypothyroidism, according 
to previous osteoporosis history. 

Methods: We conducted a cohort study from the Korean Health Insurance Review and Assessment Service claims database from Jan- 
uary 2005 to June 2006. The study population comprised women aged >65 years who had been diagnosed with hypothyroidism and 
prescribed levothyroxine monotherapy. We excluded patients who met any of the following criteria: previous fracture history, hyper- 
thyroidism, thyroid cancer, or pituitary disorder; low levothyroxine adherence; or a follow-up period <90 days. We categorized the 
daily levothyroxine doses into 4 groups: <50 ug/d, 51 to 100 ug/d, 101 to 150 ug/d, and > 1 50 ug/d. The hazard ratios (HRs)and95% 
confidence intervals (CIs) were estimated with the Cox proportional hazard model, and subgroup analyses were performed according 
to the osteoporosis history and osteoporosis-specific drug prescription status. 

Results: Among 1 1 155 cohort participants, 35.6% had previous histories of osteoporosis. The adjusted HR of fracture for the > 150 ug/ 
d group, compared with the 51 to 100 ug/d group, was 1.56 (95% CI, 1.03 to 2.37) in osteoporosis subgroup. In the highly probable 
osteoporosis subgroup, restricted to patients who were concurrently prescribed osteoporosis-specific drugs, the adjusted HR of frac- 
ture for the > 1 50 ug/d group, compared with the 51 to 100 ug/d group, was 1.93 (95% CI, 1.14 to 3.26). 

Conclusions: While further studies are needed, physicians should be concerned about potential levothyroxine overtreatment in el- 
derly osteoporosis patients. 
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INTRODUCTION 

Hypothyroidism is common in elderly women, with preva- 
lence rates of 3% to 10% in the general population and 11% 
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to 15% in the elderly population [1]. In Korea, a cohort study 
by Choi et al. [2] reported that the subclinical hypothyroidism 
prevalence was 18.9% in elderly Korean women. Therefore, 
more than 20% of elderly patients receive long-term levothy- 
roxine replacement therapy [1]. Despite the importance of this 
treatment, previous studies on the association between levo- 
thyroxine use and fracture risk have reported controversial re- 
sults [3-5]. 

During normal aging, regulatory functions such as thyroid 
hormone production and degradation decrease, and therefore 
the required doses for elderly patients with hypothyroidism 
differ from those for younger patients [6,7]. Chronic under or 
over-replacement is common in clinical practice, and data in- 
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dicate that over-replacement occurs in approximately 20% of 
levothyroxine patients who are treated with levothyroxine due 
to non-adherence or a lack of levothyroxine replacement moni- 
toring [8-10]. Bone loss due to levothyroxine over-replacement 
frequently occurs in postmenopausal women [11,12]. However, 
evidence of an association between levothyroxine doses and 
fracture risk has been insufficient. 

A recent epidemiological study by Turner et al. [1 3] reported 
that high levothyroxine doses increased the fracture risk by 
3-fold over that associated with low doses in elderly hypothy- 
roidism patients. However, the authors did not consider frac- 
ture risk stratification according to the patients' osteoporosis 
statuses. Furthermore, the American Association of Clinical En- 
docrinologists (AACE) and American Thyroid Association (ATA) 
guidelines for hypothyroidism management did not recom- 
mend contraindicated levothyroxine doses to minimize frac- 
ture risks according to the osteoporosis status [14]. 

Additionally, a study on the association between fracture 
risk and levothyroxine use in Korean patients has not yet been 
performed. Therefore, the objective of our cohort study, which 
employed a nationwide claims database, was to evaluate the 
association between levothyroxine dosage and fracture risk. 
Furthermore, we evaluated differences in this association ac- 
cording to the degree of osteoporosis. 

METHODS 

Data Source and Ethical Considerations 

We accessed the Korean Health Insurance Review and As- 
sessment Service (HIRA) database from January 2005 to June 
2006. This nationwide database contains information about 4 
159 309 elderly patients aged >65 years and 100 838 744 
prescriptions. The claims database comprises inpatient and 
outpatient care information from all clinics and hospitals, in- 
cluding demographic information, diagnoses coded according 
to the International Classification of Disease, tenth revision 
(ICD-1 0), and prescription records such as the generic and brand 
names of drugs, prescription dates, durations, dosages, costs, 
and routes of administration [15]. This database has been used 
for many other epidemiological studies [1 6-1 8]. 

From the HIRA database, we received data about elderly pa- 
tients (aged >65 years) that were extracted and de-identified 
by HIRA to protect privacy, according to the Act on the Protec- 
tion of Maintained Personal Information. The study was ex- 
empted from review by the institutional review board of the 



Seoul National University College of Medicine and Seoul Na- 
tional University Hospital because a de-identified secondary 
database was used. 

Cohort Construction 

The cohort participants were defined as elderly patients 
(aged >65 years) who were diagnosed with hypothyroidism 
(ICD-1 0: E01-E03) and prescribed thyroid hormone therapy 
from January 2005 to June 2006. Our cohort participants were 
restricted to elderly women. To select only patients who were 
prescribed levothyroxine monotherapy, we excluded patients 
who were prescribed liothyronine or liothyronine/levothyrox- 
ine combination therapy. We also excluded patients who met 
any of the following criteria: a past history of fracture (ICD-1 0: 
S22.0, S22.1, S32, S42.2-S42.4, S42.7-S42.9, S52, S62.0-S62.4, 
S62.8, S72, S82) prior to the first levothyroxine prescription 
[13]; a past history of hyperthyroidism (ICD-1 0: E05), thyroid 
cancer (ICD-1 0: C73), or pituitary disorder (ICD-1 0: E22, E23) 
[13]; low adherence to levothyroxine therapy; or a follow up 
period <90 days. Low adherence to levothyroxine therapy 
was defined by the proportion of days covered (PDC) and a 
maximum levothyroxine dose criterion of <80% and >300 
ug/d, respectively [19,20]. PDC was calculated as the number 
of days covered by prescription, divided by the number of to- 
tal follow-up days [21]. 

Thus, we selected hypothyroidism patients who adhered 
well to the prescribed levothyroxine monotherapy and had no 
previous history of any type of fracture, hyperthyroidism, thy- 
roid cancer, or pituitary disorder. The cohort participants were 
followed up from the date of first levothyroxine prescription 
within the study period to the date of the first fracture event, 
death (ICD-1 0: R96, R98, R99, 1461 ), or the last date of the study 
period (June 30, 2006), whichever occurred first. 

Exposure and Outcome Assessment 

In Korea, a total of 1 8 levothyroxine formulas were introduced 
into the market from 2005 to 2006, and these were classified 
into 2 doses, 50 and 100 ug. We calculated the mean daily le- 
vothyroxine dose as the sum of each daily levothyroxine dose 
during the follow-up period, divided by the number of total 
follow-up days. If overlap occurred between prior and later 
prescriptions, we used the later prescription dose. As the most 
commonly prescribed dose was 100 ug/d and, generally, dose 
modifications in clinical situations were 25 ug/d, we catego- 
rized the daily levothyroxine doses into 4 groups: <50 ug/d, 
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51 to 100 ug/d, 101 to 150 ug/d, and >150 ug/d. Furthermore, 
the 51 to 100 ug/d group was considered a reference category 
when hazard ratios (HRs) were calculated. 

The first fracture diagnosed in the department of orthope- 
dic surgery or emergency medicine after levothyroxine pre- 
scription was considered an outcome event. A fracture associ- 
ated with thyroid hormone use was defined as a fracture of 
the wrist or forearm (ICD-10: S52, S62.0-S62.4, S62.8), shoulder 
or upper arm (ICD-1 0: S42.2-S42.4, S42.7-S42.9), thoracic spine 
(ICD-10: S22.0, S22.1), lumbar spine or pelvis (ICD-1 0: S32), hip 
or femur (ICD-10: S72), or lower leg or ankle (ICD-10: S82) [13]. 
Patients who concurrently experienced seizure (ICD-10: G40, 
G41, R56), trauma (ICD-1 0:T79,T98,V01-V99, W11-W17, W50- 
W52, W64, X34-X39, Y01-Y04, Y30-Y32), bone malignancy (ICD- 
10: C40, C41), multiple myeloma (ICD-10: C90), or pathological 
fracture (ICD-10: M84.4, M90.7) at any time during the study 
period were excluded from the outcome definition and were 
censored at the day on which the event occurred. 

Statistical Analyses 

We evaluated the characteristics of the cohort participants, 
including age, comorbidities, Charlson comorbidity scores, co- 
medications, and the number of health service uses, with the 
chi-squared test and analysis of variance (ANOVA) as appropri- 
ate. If the descriptive analyses for each variable were shown to 
differ significantly between the groups, we performed an ad- 
ditional Cochran-Armitage trend test for categorical variables 
or Duncan's test for continuous variables. 

During the study period, we considered comorbidities and 
comedications as the confounding factors of the association 
between levothyroxine use and fracture risk. Comorbidities 
possibly related to the fracture risk were identified in claims by 
the ICD-10 codes. These comorbidities included osteoporosis 
(M81, M82), malignant neoplasm (C), diabetes mellitus (DM, 
E10-E14), hypoparathyroidism (E20), Cushing's syndrome (E24), 
adrenal disorder (E25-E27), chronic obstructive pulmonary 
disease (COPD, J40-J44), asthma (J45), chronic liver disease 
(B18, B19, K70-K77), chronic renal failure (N18, N19), heart fail- 
ure (I50), ischemic heart disease (I20-I25), arrhythmia (I44-I49), 
dementia (F00-F03), stroke (I60-I64), and epilepsy (G40, G41). 
Comedications such as corticosteroids, anti-coagulants, anti- 
epileptics, anti-depressants, benzodiazepines, proton pump 
inhibitors (PPI), and thiazolidinediones were considered po- 
tential risk factors for fracture. In contrast, bisphosphonates 
and raloxifene were considered protective factors against frac- 



ture [22]. 

The Cox proportional hazard model was used to calculate 
the HRs and 95% confidence intervals (Cls).The above-listed 
factors were considered possible confounders in the multivari- 
ate models because they were well-known fracture risk factors. 
When analyzing the fracture risk, we considered the left trun- 
cation of data, using the SAS survival data analysis function 
[23]. We performed subgroup analyses according to the previ- 
ous histories of osteoporosis diagnoses and accordingly divid- 
ed the patients into 2 subgroups: the osteoporosis subgroup 
and the non-osteoporosis subgroup. The osteoporosis sub- 
group was further divided into 2 groups according to the os- 
teoporosis-specific drug prescription: the highly probable os- 
teoporosis subgroup included patients who were prescribed 
bisphosphonate or raloxifene, while the probable osteoporosis 
subgroup included patients who were not prescribed bisphos- 
phonate or raloxifene. These operational definitions of osteo- 
porosis statuses in our study were determined from an algo- 
rithm that was formerly applied to ascertainments of osteopo- 
rosis in other HIRA database studies [1 8]. 

We calculated the incidence rates among the total elderly 
women in our dataset, the eligible cohort participants, and 
each subgroup stratified by osteoporosis statuses and daily le- 
vothyroxine dose groups. Furthermore, these incidence rates 
were standardized to the age distribution of all Korean elderly 
women in 2005. 

All statistical analyses were performed with the SAS version 
9.3 (SAS Inc., Cary, NC, USA). 

RESULTS 

Study Population 

Of the 4 159 309 elderly patients in the database, 81 542 
were diagnosed with hypothyroidism during the study period, 
63 998 of whom were women. Among the hypothyroidism 
patients, 32 743 women were prescribed thyroid hormone 
therapy. We selected the 29 081 women who were prescribed 
levothyroxine monotherapy by excluding 3662 women who 
were prescribed liothyronine monotherapy or liothyronine/le- 
vothyroxine combination therapy. Next, we applied the exclu- 
sion criteria, and 11 155 women were finally included as co- 
hort participants (Figure 1). Among the cohort participants, 
the median and interquartile ranges of follow-up duration 
were 514 days and 451 to 533 days. Among those, 256 patients 
experienced fractures during the study period. 
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4 159 309 Individuals in the HIRA database who were aged >65 years 



81 542 Individuals who were diagnosed with hypothyroidism 



63 998 Women who were diagnosed with hypothyroidism 



32 743 Women with thyroid hormone prescriptions 



3662 Women receiving LT3 monotherapy or LT3/LT4 combination therapy 



29 081 Women receiving LT4 monotherapy 



Excluded patients who met any of the following conditions: 
22 Women with a history of fracture before using levothyroxine 
7251 Women with a history of hyperthyroidism, thyroid cancer, or pituitary disorder 
2072 Women who were followed up for <90 days 
12 341 Women whose levothyroxine PDC during the study period <80% 
222 Women who received an abnormal daily levothyroxine dose range (>300 ug/d) 



Eligible participants: 11 155 women 



Figure 1. Flow chart of participant selection. HIRA, Korean Health Insurance Review and Assessment Service; LT3, liothyronine; 
LT4, levothyroxine; PDC, proportion of days covered. 



Baseline Characteristics of the Participants 
According to Daily Levothyroxine Dose and 
Osteoporosis History 

Table 1 shows the differences in the baseline characteristics 
among the daily levothyroxine dose groups. Regarding the 
daily levothyroxine doses among the total cohort participants, 
the median was 95.7 ug/d, and the interquartile range was 
81.3 to 100.0 ug/d. After categorization, 6097 patients (54.7%) 
were prescribed doses between 51 ug/d and 100 ug/d. Further- 
more, 2296 (20.6%), 1674 (15.0%), and 1088 (9.8%) partici- 
pants were prescribed daily levothyroxine doses of <50 ug/d, 
101 to 150 ug/d, and > 1 50 ug/d, respectively. The patients in 
the lowest dose group (< 50 ug/d) had a higher prevalence of 
cardiovascular problems, including heart failure, ischemic heart 
failure, arrhythmia, and stroke, than did those in the higher 
dose groups. In contrast, the highest dose group ( > 1 50 ug/d) 
had a higher prevalence of malignancy and hypoparathyroid- 
ism. However, the rates of previous histories of osteoporosis 
diagnosis and bisphosphonate prescriptions were not signifi- 
cantly different between the levothyroxine dose groups (p= 
0.46 and p = 0.36, respectively). 

We analyzed differences in the baseline characteristics of 
the subgroups according to the previous histories of osteopo- 
rosis diagnosis and osteoporosis-specific drug prescriptions 



(Table 2). Individuals in the highly probable osteoporosis sub- 
group had a higher prevalence of Cushing syndrome, adrenal 
disorders, COPD, and asthma, higher Charlson comorbidity 
scores, and more prescriptions for glucocorticoids, anti-epilep- 
tics, anti-depressants, benzodiazepines, and PPI, as well as 
higher numbers of health service uses than those in other 
subgroups, and all of the p-values for the trends associated 
with the above-listed factors were <0.01 . In contrast, the non- 
osteoporosis subgroup had a higher prevalence of chronic 
kidney diseases and heart failure, and the p-values for these 
trends were <0.01 and 0.02, respectively. 

Association Between the Daily Levothyroxine 
Dose and Fracture Risk 

Compared to the reference group (51 to 100 ug/d), the frac- 
ture risks in the <50 ug/d, 101 to 150 ug/d, and >150 ug/d 
groups were not significant in the total cohort participants. In 
a subgroup analysis according to osteoporosis diagnoses, the 
fracture risk of the highest daily dose group (> 1 50 ug/d) was 
significantly higher than that of the reference group in the os- 
teoporosis subgroup (adjusted HR, 1.56; 95% CI, 1.03 to 2.37) 
(Table 3). After we divided the osteoporosis subgroup into 2 
subgroups (highly probable osteoporosis and probable osteo- 
porosis subgroup), the fracture risk of the highest daily dose 
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Table 1. Baseline characteristics of the cohort participants who were prescribed levothyroxine monotherapy among a popula- 
tion of elderly women with hypothyroidism, according to the daily levothyroxine dose 



Characteristics 


Total (%) 


< 50 |ig/d 


Daily levothyroxine dose (%) 
51-100 ug/d 101-150 ug/d 


> 150 ug/d 


p-value 


Age (mean + SD) 


70.7±5.0 


71 .3 + 5.3 


in 1 4- R n 
/ U. / X o.u 


7094-/1 c. 
/U.Z I4.D 


70.5+4.9 


<0.01 


65-69 


5510(49.4) 


1029 (44.8) 


9007 lAQ Ql 
OUO/ (4y.oJ 


QQ/l /C9 Q\ 
oo4 \bZ.Oj 


560(51.5) 




70-74 


3213(28.8) 


676 (29.4) 


1 707 /9Q C\ 

1 161 (Zo.b) 


AOA IOQ E\ 

4y4 (zy.bj 


306 (28.1) 




75-79 


1726(15.5) 


404(17.6) 


CMC M c c\ 
343 ( 1 b.bj 


99/1 /1 9 A\ 
ZZ4 ( 1 0.4) 


153(14.1) 




>80 


706 (6.3) 


187 (8.1) 


07Q ic 0\ 
o/o (□./) 


79 IA 9\ 

/Z (4.o) 


69 (6.3) 




Comorbidity 














Osteoporosis 


3976 (35.6) 


801 (34.9) 


Z I 0/ (Ob.bJ 


CO/I 101 91 
bZ4 (0/ .0) 


384 (35.3) 


0.46 


Malignancy 


639(5.7) 


112(4.9) 


OCC /C Q\ 
ODD \D.O) 


Q9 /E E\ 

yz (b.b) 


80 (7.4) 


0.03 


Diabetes mellitus 


4364 (39.1) 


933 (40.6) 


9070 100 Q\ 

zovu (oo.yj 


CCO /9Q c\ 

bbo (oy.b) 


398 (36.6) 


0.14 


Hypoparathyroidism 


116(1.0) 


15(0.7) 


c9 /n ai 
DL (u.yj 


99 /1 /I\ 
Zo ( I .4) 


26 (2.4) 


<0.01 


Cushing syndrome 


47(0.4) 


3(0.1) 


99 /n El 
oo (U.b) 


o (U.b) 


3 (0.3) 


0.06 


Adrenal disorder 


68(0.6) 


15(1.6) 


/in /n 7\ 

4U (u. I } 


1 n in ri 

I U (U.DJ 


3 (0.3) 


0.51 


COPD 


2010(18.0) 


401 (17.5) 


1 non /1 7 q\ 
i uyu \ i /.yj 


9QC /1 7 CI 

zyb ( i / .bj 


224 (20.6) 


0.13 


Asthma 


2738 (24.6) 


577 (25.1) 


1 A~IC\ 10 A 11 

1 4/U 1 ) 


one loo 7i 

030 (ZO./J 


295 (27.1) 


0.13 


Chronic liver disease 


2098(18.8) 


420(18.3) 


1 1 Q7 M Q El 

1 1 0/ ( i y.oj 


909 tM A\ 

zyz ( i / .4j 


199(18.3) 


0.22 


Chronic kidney disease 


280 (2.5) 


78 (3.4) 


1 99 iO 91 
lOO (Z.Z) 


/in /9 /ii 
4U (Z.4J 


29 (2.7) 


0.02 


Heart failure 


868 (7.8) 


235(10.2) 


AOl tl C\\ 
4Z/ \J .UJ 


1 99 11 0\ 
1 ZZ (/ .0) 


84 (7.7) 


<0.01 


Ischemic heart disease 


2547 (22.8) 


595 (24.2) 


1 9E0 /99 0\ 
1 030 \LL.L) 


971 /99 9\ 
0/ 1 (ZZ.ZJ 


228 (21.0) 


<0.01 


Arrhythmia 


906 (8.1) 


250(11.1) 


MO (1 1\ 
4/ Z \l .1) 


mn /r ni 

I UU (D.UJ 


84 (7.7) 


<0.01 


Dementia 


325(2.9) 


95(4.1) 


1 CE 11 71 

1 bb [I. I) 


AO 10 cl 

4Z (Z.bJ 


23 (2.1) 


<0.01 


Stroke 


1256(11.3) 


297(12.9) 


COO /1 1 01 

boo ( I 1 .0] 


1 CO /O 1 1 

ibo (y. i ) 


118(10.9) 


<0.01 


Epilepsy 


239 (2.1) 


62 (2.7) 


1 99 10 ni 


09 /1 Q\ 

jz ( i .yj 


23 (2 1 ) 


0.22 


Charlson comorbidity score (mean ± SD) 


2.23 + 1.94 


2.40 + 1.93 


9 n _i_ 1 nn 
Z. I / ± I .yu 


9 iq-lo ni 
Z. I o± Z.U I 


2.29 + 2.03 


<0.01 


0 


1952(17.5) 


348(15.2) 


1 nQ9 /1 7 Q\ 

\ uyz ( i / .yj 


01 Q /1 q n\ 

o i o \ i y.uj 


194(17.8) 




1 


2722 (24.4) 


496 (21.6) 


1 ERR /9E 1\ 

\ boo \ld.i } 


A 1 9 /9/I Rl 
4 I Z \Z4.bJ 


248 (22.8) 




2 


2388 (21.4) 


508 (22.1) 


1 977 /9n Q\ 

i Li i (zu.yj 


OCO /91 1 \ 
ODO \Z 1 . 1 1 


250 (23.0) 




3 


1756(15.7) 


391 (17.0) 


0C1 MC CI 

yb I ( I b.b) 


OCO M E "71 

ZbZ (lb./) 


152(14.0) 




>4 


2337 (21.0) 


553 (24.1) 


1911 HQ CI 

I Z I I (I y.b) 


090 (10 71 

jzy ( i y./j 


244 (22.4) 




Co-medication 














Glucocorticoid 


6614(59.3) 


1344(58.5) 


OCC1 /CO 01 

obbl (bo.Z) 


1 noo /ci /i\ 
I UZo (b 1 .4) 


691 (63.5) 


<0.01 


Anti-coagulants 


266 (2.4) 


76 (3.3) 


1 9Q (O 11 
1 Zo (Z. \ } 


/in /9 a\ 

4U (Z.4J 


22 (2.0) 


0.01 


Anti-epileptics 


1520 (13.6) 


325(14.2) 


□no / 1 9 91 
oUO ( 1 O.ZJ 


91 c M 9 Q\ 
Z I b ( I Z.O) 


177 (16.3) 


0.03 


Anti-depressants 


2446 (21.9) 


538 (23.4) 


191/1 /91 c\ 
I 0 1 4 (Z I 


0/I9 /9n /I\ 
04Z (ZU.4J 


252 (23.2) 


0.08 


Benzodiazepines 


7327 (65.7) 


1521 (66.3) 


9009 /RE E\ 

jyyz \bb.bj 


1 non /rc 1 \ 
I uyu (bb. I J 


724 (66.5) 


0.79 


Proton pump inhibitor 


1153(10.3) 


266(11.6) 


RfiZL fQ Ql 

UU4- \CJ.CJ/ 


I Jj W.J/ 


124(11.4) 


0.05 


Thiazolidinedion 


371 (3.3) 


89 (3.9) 


190 (3.1) 


54 (3.2) 


38 (3.5) 


0.37 


Raloxifene 


351 (3.2) 


76 (3.3) 


202 (3.3) 


50 (3.0) 


23 (2.1) 


0.19 


Bisphosphonate 


2089(18.7) 


436(19.0) 


1151 (18.9) 


320(19.1) 


182(16.7) 


0.36 


No. of health service use (mean + SD) 


34.0 + 17.0 


34.3 + 16.9 


33.7 + 16.9 


33.7 + 16.9 


35.2+18.0 


<0.01 


1-19 


2119(19.0) 


415(18.1) 


1167(19.1) 


325(19.4) 


212(19.5) 




20-39 


5564 (49.9) 


1156(50.4) 


3066 (50.3) 


835 (49.9) 


507 (46.6) 




40-59 


2585 (23.2) 


550 (24.0) 


1392 (22.8) 


379 (22.6) 


264 (24.3) 




>60 


887(8.0) 


175(7.6) 


472 (7.7) 


135(8.1) 


105(9.7) 




Total 


11 155(100.0) 


2296(100.0) 


6097(100.0) 


1674(100.0) 


1088(100.0) 





COPD, chronic obstructive pulmonary disease. 
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Table 2. Baseline characteristics of subgroup participants who were prescribed levothyroxine monotherapy among a population 
of elderly women with hypothyroidism, according to the osteoporosis status 



Characteristics 


Highly probable osteoporosis (%)' 


Probable osteoporosis (%) 2 


Non-osteoporosis (%) 3 


p-value 


Daily levothyroxine dose (mean ± SD, pg/d) 


98.2+42.2 


101.6 + 44.5 


99.0+43.5 


0.02 


<50 


453 (21.0) 


468(19.2) 


1375(21.0) 




51-100 


1193(55.3) 


1306(53.6) 


3598 (54.8) 




101-150 


328(15.2) 


404(16.6) 


942(14.4) 




>150 


183 (8.5) 


259(10.6) 


646 (9.9) 




Age (mean ±SD) 


71.0 + 4.9 


70.1+4.6 


70.8 + 5.2 


<0.01 


65-69 


987 (45.8) 


1306(53.6) 


3217(49.0) 




70-74 


676 (31.3) 


699 (28.7) 


1838 (28.0) 




75-79 


360(16.7) 


331 (13.6) 


1035(15.8) 




>80 


134 (6.2) 


101 (4.1) 


471 (7.2) 




Comorbidity 










Malignancy 


140 (6.5) 


138(5.7) 


361 (5.5) 


0.23 


Diabetes Mellitus 


829 (38.4) 


965 (39.6) 


2570 (39.2) 


0.72 


Hypoparathyroidism 


24(1.1) 


60 (2.5) 


32 (0.5) 


<0.01 


Cushing syndrome 


20 (0.9) 


7(0.3) 


20 (0.3) 


<0.01 


Adrenal disorder 


30(1.4) 


10(0.4) 


28 (0.4) 


<0.01 


COPD 


444 (20.6) 


479(19.7) 


1087(16.6) 


<0.01 


Asthma 


591 (27.4) 


659 (27.0) 


1488 (22.7) 


<0.01 


Chronic liver disease 


402(18.6) 


527 (21.6) 


1169(17.8) 


<0.01 


Chronic kidney disease 


32(1.5) 


41 (1.7) 


207 (3.2) 


<0.01 


Heart failure 


150 (7.0) 


168 (6.9) 


550 (8.4) 


0.02 


Ischemic heart disease 


508 (23.6) 


566 (23.3) 


1473 (22.5) 


0.50 


Arrhythmia 


174 (8.1) 


175(7.2) 


557 (8.5) 


0.13 


Dementia 


77 (3.6) 


64 (2.6) 


184 (2.8) 


0.12 


Stroke 


263(12.2) 


277(11.4) 


716(10.9) 


0.26 


Epilepsy 


54 (2.5) 


51 (2.1) 


134 (2.0) 


0.43 


Charlson comorbidity score (mean + SD) 


2.41 ±1.95 


2.31+1.91 


2.14 + 1.94 


<0.01 


0 


333(15.4) 


383(15.7) 


1236(18.8) 




1 


454 (21.1) 


556 (22.8) 


1712(26.1) 




2 


462 (21.4) 


550 (22.6) 


1376 (21.0) 




3 


398(18.5) 


413(17.0) 


945(14.4) 




>4 


510(23.6) 


535 (22.0) 


1292(19.7) 




Co-medication 










Glucocorticoid 


1462 (67.8) 


1593(65.4) 


3559 (54.2) 


<0.01 


Anti-coagulants 


47 (2.2) 


49 (2.0) 


170 (2.6) 


0.22 


Anti-epileptics 


408(18.9) 


347(14.2) 


7665(11.7) 


<0.01 


Anti-depressants 


546 (25.3) 


563 (23.1) 


1337 (20.4) 


<0.01 


Benzodiazepines 


1538 (71.3) 


1675 (68.7) 


4114(62.7) 


<0.01 


Proton pump inhibitor 


279(12.9) 


299(12.3) 


575 (8.8) 


<0.01 


Thiazolidinedion 


69 (3.2) 


73 (3.0) 


229 (3.5) 


0.48 


No of health service use (mean ± SD) 


38.8+18.7 


37.8 + 17.8 


30.9 + 15.4 


<0.01 


1-19 


270(12.5) 


318(13.1) 


1531 (23.3) 




20-39 


965 (44.7) 


1150 (47.2) 


3449 (52.6) 




40-59 


640 (29.7) 


680 (27.9) 


1265(19.3) 




>60 


282(13.1) 


289(11.9) 


316(4.8) 




Total 


2157(100.0) 


2437(100.0) 


6561 (100.0) 





COPD, chronic obstructive pulmonary disease. 

'Included patients who were diagnosed with osteoporosis and were prescribed bisphosphonate or raloxifene, 
included patients who were diagnosed with osteoporosis, but were not prescribed bisphosphonate or raloxifene. 
'Included patients who were neither diagnosed with osteoporosis nor prescribed bisphosphonate. 
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Table 3. Association between the daily levothyroxine dose and the fracture risk, according to the osteoporosis status in elderly 
women with hypothyroidism 



Daily levothyroxine 
dose (|ig/d) 


No. of participants (%) 


No. of events 


Incidence rate 
</1 000 PY) 


cHR (95% CI) 


aHR 1 (95% CI) 


Osteoporosis 2 












<50 


921 (20.1) 


4b 


42.4 


1 .11 (U.ob, 1 ./b) 


1 1 c /n on 1 cc\ 
1. lb (U.oU, l.bb) 


51-100 


2499 (54.4) 


1 1 c 
lb 


ob.U 


Reference 


Reference 


101-150 


732(15.9) 


37 


37.9 


1.09 (0.75, 1.58) 


1.16(0.80, 1.69) 


>150 


442 (9.6) 
4594(100.0) 


7Q 

111 


cn 7 

OJ.U 


1 AQ (D QQ 7 1A\ 
I .HU \U.DZJ, Z.Z4/ 


1 .00 ( 1 .Do, Z. J / / 


Highly probable osteoporosis 3 










<50 


453 (21 0) 


LL 


A 1 □ 


U.y/ (U.bU, I .□/) 


n oi tr\ cc 1 r-rn 
U.30 (U.bb, I .bb) 


51-100 


1193(55.3) 


0/ 


A A 1 

44. o 


Reference 


Reference 


101-150 


328(15.2) 


18 


41.1 


0.91 (0.54, 1.54) 


0.92 (0.54, 1.58) 


>150 


183 (8.5) 
2157(100.0) 


1 Q 

i a 

1 9K. 
I ZD 


P.1 Q 

o I .a 

/IP. C 
40. u 


i oc M 1 1 q m\ 

1 .00 l 1 . 1 1 , J.U/ J 


i qo|i 14 i ic\ 

1 ,3J \ 1 . 1 H, J. ZD/ 


Probable osteoporosis 4 












<50 


468(19.2) 


l/l 

Z4 


ai n 


1 co /n m o aq\ 
I .bZ (U.ldo, Z.4o) 


1 AO /fl O/l T /in 

1 .4o (U.o4, Z.4Z) 


51-100 


1306(53.6) 


AQ 
4o 


1Q E 
Zo.D 


Reference 


Reference 


101-150 


404(16.6) 


19 


35.3 


1.23(0.72, 2.09) 


1.60 (0.90, 2.82) 


>150 


259(10.6) 
2437(100.0) 


10 
101 


29.4 
32.4 


1 03 10 2 03) 




Non-osteoporosis 5 












<50 


1375(21.0) 


52 


31.7 


1.07(0.78, 1.47) 


1.05(0.75, 1.47) 


51-100 


3598 (54.8) 


138 


29.9 


Reference 


Reference 


101-150 


942(14.4) 


44 


35.4 


1.18(0.84, 1.66) 


1.10(0.78, 1.56) 


>150 


646 (9.9) 
6561 (100.0) 


22 
256 


26.2 
30.7 


0.88(0.56, 1.37) 


0.84(0.52, 1.35) 



2 533 289 



175 783 



46.3 



PY, person-year; cHR, crude hazard ratio; CI, confidence intervals; aHR, adjusted hazard ratio. 
'Adjusted for age, comorbidities, comedications, Charlson comorbidity score, and number of health service uses, 
included patients who were diagnosed with osteoporosis, regardless of bisphosphonate or raloxifene prescription status. 
'Included patients who were diagnosed with osteoporosis and were prescribed bisphosphonate or raloxifene. 
'Included patients who were diagnosed with osteoporosis, but were not prescribed bisphosphonate or raloxifene, 
included patients who were neither diagnosed with osteoporosis nor prescribed bisphosphonate. 



group was significantly increased only in the highly probable 
osteoporosis subgroup (adjusted HR, 1.93; 95% CI, 1.1 4 to 3.26). 
Furthermore, the fracture risk trend in the highly probable os- 
teoporosis subgroup was correlated linearly with the daily le- 
vothyroxine dose (p for trend: 0.03). 

DISCUSSION 

In this nationwide retrospective cohort study on the associa- 
tion between the levothyroxine dose and fracture risk among 
elderly women (aged >65 years), we found a significant asso- 
ciation between a higher dose of levothyroxine (> 1 50 ug/d) 
and fracture risk in the highly probable osteoporosis subgroup. 



This result suggests that high-dose levothyroxine treatment 
might increase the fracture risk in severely osteoporotic pa- 
tients. However, a levothyroxine treatment < 1 50 ug/d was 
not associated with a fracture risk, regardless of osteoporosis 
status. 

When we compared the age-standardized incidence rates 
among levothyroxine doses of > 1 50 ug/d in the highly proba- 
ble osteoporosis subgroup, eligible cohort participants, and all 
elderly women in our dataset, the incidence rate ratios of the 
> 1 50 ug/d levothyroxine group were 2.54 when compared to 
the eligible cohort participants and 1 .72 when compared to all 
elderly women in our dataset. 

The significantly increased fracture risk in severely osteopo- 
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rotic patients agreed with those reported in previous studies. 
Evidence about the association between levothyroxine replace- 
ment and reduced bone mineral density (BMD) was well sup- 
ported. Faberand Galloe [11] reported in a meta-analysis of 13 
studies that women with an average age of 39.6 years who 
were treated with 1 64 ug of levothyroxine per day had a 2.67% 
lower BMD than that of the controls. Uzzan et al. [12] also re- 
ported in a meta-analysis of 25 studies that bone losses 
caused by long-term levothyroxine use in postmenopausal 
women significantly decreased BMDs by 7% in the lumbar 
spine and 9% in the femoral neck (mean age, 61.1 years; mean 
follow-up duration, 9.6 years). Furthermore, several research- 
ers reported that patients who received levothyroxine replace- 
ment therapy, especially postmenopausal women, had an in- 
creased fracture risk. Flynn et al. [4] conducted a retrospective 
cohort study of the association between the occurrence of ad- 
verse outcomes such as cardiac and osteoporotic events and 
thyroid-stimulating hormone (TSH) concentrations in patients 
who were prescribed levothyroxine replacement therapy. The 
TSH-suppressed patients, represented as high-dose levothyrox- 
ine prescriptions, had a 2.02-fold increase (95% CI, 1.55 to 2.62) 
in the risk of osteoporotic fractures. In a nested case-control 
study conducted in Canada by Turner et al. [13], 213 511 levo- 
thyroxine users aged 70 years or older from a health insurance 
database were followed up for 3.8 years, and the odds ratio for 
fracture was 3.45 (95% CI, 3.27 to 3.65) for the highest dose 
group (>93 ug/d), compared to the lowest dose group (<44 
ug/d). This study showed that current levothyroxine treatment 
was associated with a significantly increased fracture risk in a 
strong dose-response relationship. Since previous studies did 
not consider the patients' baseline osteoporosis statuses, fur- 
ther studies are needed to evaluate the safe levothyroxine 
ranges according to the osteoporosis status. 

These results agreed with the pathologic mechanisms of iat- 
rogenic hyperthyroidism caused by excessive levothyroxine 
use [24]. The pathogenic mechanism that affects the bones in 
hyperthyroidism is based on an increase in both the number 
and turnover rate of bone turnover units and thus increases in 
osteoclast and osteoblast activity, with 50% reductions in the 
remodeling cycle and increased unit activation frequency. 
These changes lead to an uncoupling of resorption and forma- 
tion, and the net result of mineralized bone losses in varying 
amounts depends on factors such as sex, menstrual status, 
thyroid disease severity, and the sum of other osteoporosis 
risk factors [25,26]. The result of a previous study about the in- 



teraction between iatrogenic hyperthyroidism-mediated re- 
duced BMD and osteoporosis severity were limited. However, 
these 2 conditions should interact strongly with each other 
because they share the similar pathologic mechanism of re- 
duced BMD. Additionally, in our study, an significant interac- 
tion between the levothyroxine dose and osteoporosis status 
was observed in logistic regression (p<0.001). 

Our results indicated no significant relationship between 
the daily levothyroxine dose and fracture risk in non-severe 
osteoporosis patients. When we calculated the appropriate 
sample size assuming 80% power, 5% alpha-error, and a 3% 
expected fracture incidence in unexposed group, 1492 patients 
were required to confirm a relative risk of 2.0. In our study, the 
minimum number of participants in 1 comparison was 1565 
patients in the non-severe osteoporosis subgroup. Therefore, 
the number of participants in this study's subgroup analyses 
was sufficient to examine the relationship between the daily 
levothyroxine dose and fracture risk. 

We considered the potential occurrence of "confounding by 
indication," in which prognostic factors influence the treatment 
decisions. In our study, we evaluated the factors associated 
with the indications used to decide the levothyroxine doses 
and fracture risk factors. The AACE and ATA guidelines for hy- 
pothyroidism management suggested that patients with an- 
gina symptoms should receive half of the regular levothyrox- 
ine dose at the start of levothyroxine treatment for hypothy- 
roidism; furthermore, the target TSH level was also higher in 
this population [14]. To evaluate the previous medical histories 
associated with angina symptoms, we analyzed covariates 
such as ischemic heart disease, arrhythmia, stroke, and heart 
failure in multiple logistic regressions. However, we did not 
find any associations between the fracture risk and these vari- 
ables in any osteoporosis subgroup. Furthermore, we evaluat- 
ed the differences in the average daily levothyroxine doses in 
the lowest dose groups according to the degree of osteoporo- 
sis, and no significant differences were observed in our data 
with ANOVA (p=0.43; average daily levothyroxine doses: 46.5 
ug/d for the non-osteoporosis subgroup, 46.7 ug/d for the 
probable osteoporosis subgroup, and 47.2 ug/d for the highly 
probable osteoporosis subgroup). 

Another considerable bias was misclassification. According 
to the previous osteoporosis history, it was expected that pa- 
tients who were included in the highly probable osteoporosis 
subgroup would have homogenous characteristics. Until 201 1, 
physicians in Korea could prescribe osteoporosis-specific drugs 
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such as bisphosphonate and raloxifene when theT-scores of 
the lumbar spine or femur neck were less than -3 according to 
dual-energy X-ray absorptiometry (DXA) because of unitary 
insurance coverage regulations. Therefore, a bisphosphonate 
or raloxifene prescription indicated that a patient's T-scores ac- 
cording to DXA were less than -3, consequently, the highly 
probable osteoporosis subgroup comprised homogenous pa- 
tients with severe osteoporosis [27]. In many epidemiologic 
studies that used claims databases, operational definitions of 
osteoporosis diagnosis such as combined ICD-10 codes of os- 
teoporosis and the application of exclusive medications for os- 
teoporosis treatment were used [1 6,28] 

In clinical practice, low patient adherence is a common 
problem associated with levothyroxine therapy [29]. To elimi- 
nate the effects of low patient adherence, we included patients 
with adequate levothyroxine adherence, according to the PDC 
(>80%) and the maximum daily dose criterion (<300 ug/d). 
PDC was developed to evaluate prescription adherence in 
chronic drug use situations [21]. Generally, PDC above 80% is 
considered appropriate for drugs used to treat chronic diseas- 
es such as hypertension and DM [20]. Additionally, we restrict- 
ed our analytical cohort to patients who were prescribed levo- 
thyroxine doses <300 ug/d. Garber et al. [14] and Ain et al. 
[1 9] reported that most patients who were regularly prescribed 
levothyroxine doses >300 ug/d had low levothyroxine adher- 
ence. 

Our results should be considered in light of several limita- 
tions. First, the data in our study did not include body weight 
and lifestyle factors such as smoking history and alcohol use 
status. However, in previous epidemiologic studies, the frac- 
ture risk estimates were little affected by adjustments for body 
weight or lifestyle factors [30,31]. Second, our study partici- 
pants included previous levothyroxine users and thus did not 
comprise wholly new users. In a sensitivity analysis restricted 
to new users (defined as not having used levothyroxine during 
the first 6 months of the study period), we identified the frac- 
ture risk according to the daily levothyroxine dose. These new 
users included 1396 patients, among whom 21 patients expe- 
rienced fractures. The HRs among the new users were 1 .89 (95% 
CI, 0.66 to 5.40), 0.58 (95% CI, 0.08 to 4.43), and 2.10 (95% CI, 
0.28 to 15.61) in the <50 ug/d, 101 to 150 ug/d, and >150 
ug/d groups, respectively, when compared to the 51 to 100 
ug/d group. Given the lack of statistical power, we could not 
conclude the fracture risk according to the levothyroxine dose 
among new users. Therefore, further studies based on longer 



follow-up periods are recommended to confirm the associa- 
tion between the levothyroxine dose and fracture risk. Third, 
the data in our study were left-truncated because the dataset 
did not contain the medical records of participants before 2005. 
Consequently, we could not ascertain the previous levothyrox- 
ine prescription and medical history statuses before the start- 
ing point of the data. Therefore, we applied the statistical meth- 
od for left-truncated data reported by Cain et al. [23]. Bias can 
be reduced when analyses account for left truncation, although 
the results are unstable when the truncated fraction is high. 
Additionally, this method adjusts for a just distribution of out- 
come events, but does not adjust for exposure. Fourth, our data 
did not contain records for the pharmaceutical drug prepara- 
tions or the actual drug ingestion of the patients in real-life. 
Therefore, we could not identify the true levothyroxine intake 
amounts of the patients. The last limitation of this study was 
related to the validity of diagnoses such as fractures and hypo- 
thyroidism. Our data did not include laboratory and clinical in- 
formation such as the BMD and TSH levels. To increase the va- 
lidity, we developed operational definitions of fracture, hypo- 
thyroidism, and osteoporosis, which were restricted to frac- 
tures diagnosed by the department of orthopedic surgery or 
emergency medicine, the diagnosis of hypothyroidism with 
thyroid hormone replacement therapy use, and the diagnosis 
of osteoporosis with osteoporosis-specific drug prescriptions, 
respectively. Additionally, a validation study was used to com- 
pare the diagnoses derived from the medical insurance claims 
database in Korea with the actual diagnoses in the medical re- 
cords. The overall positive predictive value of the diagnoses 
was 81 .8% in the overall hospitalized patient population [32]. 

Despite these limitations, our study has some strengths. The 
HIRA database was used as a data source, and accurate verifi- 
cation of the starting point, duration, and dosage of the drugs 
was available during the study period. Consequently, informa- 
tion bias, especially recall bias, was possibly minimal. Further- 
more, our novel finding provides evidence for an association 
between the fracture risk and levothyroxine dosage, according 
to the osteoporosis status. 

In conclusion, prescribed levothyroxine doses of > 150 ug/d 
were associated with a higher risk of fracture in elderly female 
patients with severe osteoporosis. While further studies are 
needed to establish evidence for the risk-benefit balance, phy- 
sicians should be concerned about potential levothyroxine 
overtreatment in elderly women with histories of severe os- 
teoporosis. 
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